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LATERAL MIGRATION AND ROTATIONAL MOTION OF ELLIPTIC PARTICLES IN 

PLANAR POISEUILLE FLOW 

DEWEI QI*, LI-SHI LUO 1 , RAJA ARAVAMUTHAN* , AND WILLIAM STRIEDER? 

Abstract. Simulations of elliptic particulate suspensions in the planar Poiseuille flow are performed by 
using the lattice Boltzmann equation. Effects of the multi-particle interaction on the lateral migration and 
rotational motion of both neutrally and non-neutrally buoyant elliptic particles are investigated. Low and 
intermediate total particle volume fraction f a = 13%, 25%, and 40% are considered in this work. 

Key words, elliptic particulate suspensions in the planar Poiseuille flow, lateral migration, rotational 
motion, multi-particle interaction, the Segre and Silberberg effect, lattice Boltzmann simulation 

Subject classification. Fluid Mechanics 

1. Introduction. Particle suspensions in pressure driven flows can be found in many industrial ap- 
plications such as transport and refining petroleum, paper manufacturing, pharmaceutical processing and 
environmental waste treatment. An understanding of the detailed dynamical behavior of particulate sus- 
pensions in pressure driven flows (lateral migration, rotational motion, and spatial distribution of particles) 
would aid designing and developing feasible and economic industrial processes. For many years the research 
activity in this area has been focused on the investigation of the dynamic behavior of suspensions in the 
Poiseuille flow. 

Segre and Silberberg [25, 26] found that a neutrally buoyant; particle migrates to an equilibrium position 
between a channel wall and the channel centerline due to the wall effect, the velocity profile curvature and 
shear force. This phenomenon, attributed to the nonlinear inertia effect, has been confirmed in both theory 
arid experiments [3, 4, 16, 8], and more recently in direct numerical simulations [6, 7, 10, 11]. 

Feng, Hu and .Joseph [6, 7] have simulated the motion of a single circular particle in the planar Poiseuille 
flow by using a finite element method. Both neutrally and non-neutrally buoyant particles were examined. 
For a neutrally buoyant particle, the Segre-Silberberg phenomenon was reproduced in the simulations and 
it was found that the equilibrium position is closer to the wall at a higher flow velocity. For a non-neutrally 
buoyant particle, when the density difference between the particle and the fluid is small, the equilibrium 
position is either close to the wall or to the centerline, depending on whether the particle leads ahead or lags 
behind the fluid velocity locally. When the density difference (or the buoyancy effect) is large enough, the 
particle always moves to the centerline regardless whether it is lighter or heavier than the fluid. In general, 
a difference in the relative velocity across a solid particle may drive the particle to move laterally since the 
side with a higher relative velocity may experience a lower pressure. Therefore, these authors suggested 
four mechanisms responsible for the particle motion in the Poiseuille flow: (1) wall lubrication repulsion, 
(2) inertial lift due to shear slip, (3) a lift due to particle rotation, and (4) a lift associated with the curvature 
of the undisturbed velocity profile. 
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Little information on the motion of elliptic particle in a pressure driven flow is available in literature. In 
the present paper, we report the direct numerical simulation results for the systems of a single elliptic particle 
or multiple elliptic particles in the planar Poiseuille flow by using the lattice Boltzmann (LB) method. Both 
neutrally and non-neutrally buoyant particles are investigated. Only low and intermediate values of the total 
particle volume (area) fraction, i.e. f a = 13%, 25%, and 40% are considered in the simulations. Effects of 
the multi-particle interaction on the migration and rotational motion of particles are also analyzed. 

2. Theory and simulation method. The fluid flow is governed by the Navier Stokes equations: 

pfdfVL + pjuVu = — VP + r/V" u + pfg (2.1) 


where pf, u, P, r), and g are the fluid density, the velocity field, the pressure field, the fluid viscosity, and 
the acceleration due to gravity, respectively. The equations of motion for the solid particle are 

dU ^ , 

in— — = mg + F, (2.2a) 

I ~ = T. (2.2b) 

dt y ’ 

where m, U, O, I, are the mass, the linear velocity, the angular velocity, and the moment of inertia tensor 
of the particle, respectively; F and T are the hydrodynamic force and torque exerted on the particle, 
respectively. 

The Navier-Stokes equations can be simulated by the lattice Boltzmann equation (LBE) [27, 5, 18, 
17, 9, 14, 15, 13, 1, 2]. The motion of nonspherical particles can be handled by solving the equations of 
motion for each particle [19, 20, 21]. In the past, the LBE method has been successfully used to simulate 
two-dimensional (2D) rectangular particles and three-dimensional (3D) beds of cylinders in sedimentation 
flows [22, 23]. In these simulations, the numerical results agree with experimental data very well and the 
direct numerical simulations by using the LBE method have quantitatively captured the essential physics of 
the particle-fluid-flow, and thus can be able to correctly predict the flow phenomena. Therefore, the same 
method is used in this work. The details of the LBE method has been reported elsewhere (cf. [14, 15, 2, 21]) 
and will not be repeated here. In the present work, it is assumed that when two particles collide, an elastic 
collision occurs. In other words, the linear and angular momenta of each particle are conserved during the 
collision. This hard shell scheme, developed by Rebertus and Sando [24] for molecular dynamic simulations 
and adopted by Qi [22, 23] for sedimenting cylinders, will be used in this work. 


3. Simulations and Results. One can find that the Poiseuille flow is determined by four dimensionless 
parameters [12] 
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where p s is the ratio of solid density p p to fluid density pf, W is the width of the channel, a is the radius of 
the principal axis of the elliptic particle, C is the confinement ratio, 7 is the local shear rate of flow, and Re 
is the particle Reynolds number. In the planar Poiseuille flow, the shear rate 
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is not a constant and varies with coordinate x in the cross flow direction, where U m is the maximum velocity 
of the undisturbed flow in the channel. We use the shear rate at x = IP/4 to estimate the particle Reynolds 
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number unless otherwise stated in this work, i.e., 


Re — 
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where v is the kinematic viscosity of the fluid. The dimensionless parameter 


S = 


Re 

~G 


2a'y 2 



(3.2) 


measures the ratio of life to buoyant weight, which will be used later. 

In the simulation presented here, a constant body force simulating the gravity is used to achieve the 
effect of constant pressure gradient in the Poiseuille flow. The flow and gravity are along the {/-direction, as 
shown in Fig. 1. And the periodic boundary conditions are applied in the {/-direction, as shown in Fig. 1. 
The bounce-back boundary conditions are applied at walls as well as the surface of particles to mimic the 
noslip boundary conditions. The bounce-back scheme used here for the fluid-solid interface of a moving 
particle is first proposed by Ladd [14], improved by Aidun el al. [2] , and is discussed in detail in [21]. 



Fig. 1. The planar Poiseuille flow configuration. Gravity and flow are in the y -direction (down). 


3.1. Neutrally buoyant elliptic single particle. Before investigating a multi-particle system, the 
migration of a singe elliptic particle is studied first. A computational domain of size 251 x 251 in lattice unit 
is used. The major radius of the ellipse is a = 15 and the minor radius is b = 7.5. The confinement ratio 
C = 0.12 for all the cases. The bulk Reynolds number of suspension is defined by 


*■ ".6 — , 

v 

where is the average velocity of the suspension and v is the kinematic viscosity. 

First, four simulations of a single neutrally buoyant elliptic particle with different initial positions in the 
^-direction (the cross flow direction) are performed at Re?, - 140.625. In each of the four runs, the elliptic 
particle migrates to the equilibrium position approximately 0.156FF away from a wall. This equilibrium 
position is independent of the initial orientation and position, as shown in Fig. 2. However, it takes a much 
longer time for a particle with an initial position close to the center of channel to reach the equilibrium 
position because the shear rate at the channel center is zero. The Segre-Silberberg phenomenon is observed 
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in the system of a neutrally buoyant elliptic particle in the planar Poiseuille flow. The results show that the 
particle always lags behind the flow locally and the effect of velocity curvature indeed pushes the particle 
toward the wall. This mechanism has been clearly explained by Feng, Hu, and Joseph [7]. 



Fig. 2. The migration of a neutrally buoyant elliptic particle with different initial positions xq. There are four simulations 
for single ellipse at Ret, — 140.625. All distances are normalized by the width of channel W . The same is applied to other 
figures. The wiggles are caused by the discrete effect of the solid particle. 

The ellipse always rotates at its equilibrium state because the particle still experiences a shear force. 
The angular velocity is periodic in time. The rotational rate at the final equilibrium state is plotted in Fig. 3 
for the neutrally buoyant particle. 



Fig. 3. The angular velocity as a Junction of time for a neutrally buoyant particle. 

In order to describe the orientation of the particle, an averaged angular distribution across the channel 
is used: 


f(6) = (6(6 - 0(t))) t > (3.3) 

where S is a Dirac function, and {•)* is the average performed over an ensemble of trajectories 6(t) at a 
given position x. The averaged angular distribution f(9) gives the probability of ellipse along angle 9. The 
angular distribution function is normalized to unity over the interval of [0, w]. Fig. 4 shows that the ellipse is 
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much more likely to align its long axis with the flow direction (y-direction) because the particle experiences 
a weaker torque from the flow when it is in alignment with the flow. 



Fig. 4. The angular distribution function of neutrally buoyant particle. 


3.2. Non-neutrally buoyant single particle. When the particle density is slightly larger than fluid 
density (p s = 1.005, Re = 6.08, and G = 27.7), the particle leads the flow velocity locally and the relative 
velocity at the side of the particle closer to the wall is increased, resulting in a lower pressure due to a higher 
relative velocity on this side. The particle migrates to a wall due to the pressure difference and stabilizes at 
a position closer to a wall than does a neutrally buoyant particle. In this case the ratio of the lift to buoyant 
weight S — 0.22. 

When the solid density is sufficiently larger than the fluid density (p s = 1.01, Re = 6.08, and G = 55.4), 
the particle motion resembles that of sedimentation. The wall repulsion overpowers the inertial lift and the 
velocity curvature effect, and the particle moves to the centerline. Figure 5 shows the migrations for the 
cases of p s = 1.005 and p s = 1.01. For the purpose of comparison, the migration for a neutrally buoyant 
ellipse (p s — 1) is also shown in the same figure. It is clear that sedimentation is the major force driving 
the particles from the wall to the channel center when the density of ellipse is large enough (e.g., p, = 1.01). 
This can be seen more clearly from the averaged angular distribution across the channel. In a sedimentation 
flow, the orientation of an ellipse along the cross flow direction is dominant due to the wakes associated with 
the nonlinear inertia effect. The angle 9 (cf. Fig. 1) of the ellipse as a function of time is depicted in Fig. 6 
for the cases of p s = 1.005 and p s = 1.01. The ellipse of p s = 1.01 turns quickly to the horizontal direction or 
the cross flow direction around 0 = 0°, while the ellipse of p s = 1.005 turns to the flow direction or vertical 
direction, oscillating approximately around 9 = 83°. It is evident that when p s = 1.01 the sedimentation 
effect is dominant since the long end of the particle turns to the horizontal direction and the particle moves 
to the centerline (x/W = 1/2), where the shear rate vanishes. In contrast, the particle with p s = 1.005 turns 
to the vertical direction and moves toward to the wall. This behavior may be attributed to the combined 
effects of the Poiseuille flow' and sedimentation. The shear lift and the wall effect force the particle to orient 
in the vertical direction. Huang, Hu and Joseph [11] pointed out that the particle may turn to the vertical 
direction due to the wall effect at a small particle Reynolds number in a sedimentation flow. Similarly, the 
wall effect also contributes to the determination of the particle orientation when p s = 1.005. In these two 
cases, the particle ceases to rotate due to the sedimentation and the wall effect. However, the rotational 
behavior of a neutrally buoyant ellipse is entirely different from that of the above two buoyant cases of p, > 1. 
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Fig. 5. The migration of ellipse with ps — 1.0, 1.005 and 1.01 at Re*, — 140.625. The wiggles are due to the discrete effect 
in the particle-fluid interfaces . 
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Fig. 6. The rotational angles of ellipse as a function of time for the cases of p s — 1.005 and 1.01 at Re& — 140.625. 


The neutrally buoyant particle is always rotating, and is much more likely to align its long axis with the 
flow direction than with the cross flow direction, as previously discussed and shown in Fig. 4. 

As previously mentioned, the neutrally buoyant particle lags slightly behind the local velocity of the 
undisturbed flow and the particle moves toward the wall. When the particle is slightly lighter than the fluid 
(p s = 0.9995, Re = 6.08, G = —2.77), the lagging velocity increases. The relative velocity at the end of the 
particle closer to a wall (cf. Fig. 1) becomes larger, and the inertial lift pushes the particle to the centerline. 
At an even lighter particle density (p s = 0.995, Re = 6.08, G — —27.7), the lagging velocity further increases, 
and a stronger wall repulsion forces the particle to move across the centerline, where the particle experiences 
the same repulsive forces in opposite direction by symmetry of the flow. Therefore, the particle swings 
back arid forth, oscillating about the centerline as shown in Fig. 7. Again, when the sedimentation effect is 
dominant for the particle of p s = 0.995, the particle ceases to rotate on average in time. Whereas for the 
particle with less buoyancy (p s = 0.9995), it never stops rotating, because the sedimentation effect is not 
dominant in this case. 
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Fig. 8. The snap shots of the particle configurations of 23 neutrally buoyant elliptic particles with the volume fraction 
fa = 13% at time steps (from left to right) t = 0, 300,000, and 100,000. 

3.3, Neutrally buoyant multi-particles. To study the multi-particle interaction, the motion of 23 
neutrally buoyant elliptic particles in the Poiseuille flow is simulated with the same computational domain. 
In this case, the particle volume (area) fraction f a is 13%. and the bulk Reynolds number is 140.625. Three 
snap shots of the particle configurations at t = 0, 300,000, and 400,000 are shown in Fig. 8. It is clearly seen 
that the particles move to the positions between the walls and the centerline. There are two maxima in the 
averaged volume fraction distribution across the channel and these maxima are located between the walls 
and centerline, as shown in Fig. 9. However, the maxima are closer to the centerline for the multi- particle 
system than for the single particle system because multi-particle interaction reduces the curvature effect of 
the flow. It seems that the particles form clusters, as shown in Fig. 8. With a careful observation, it is found 
that there are two small shoulders at the averaged velocity distribution across the channel shown in Fig. 10. 
The locations of the shoulders on the velocity distribution correspond to the maxima of the volume fraction 
distribution across the channel, indicating that the averaged velocity of particles is suppressed at areas of 
higher particle volume fraction f a . This is expected because higher effective viscosity due to higher particle 
volume fraction reduces the velocity gradient in these areas. 

Next, the averaged multi-particle angular distribution function is used to describe the orientation of 

7 




0.6 


total particle volume fraction=1 3% 



Fig. 9. The particle volume fraction distribution across the channel for neutrally buoyant particles (p s — 1.0) at f a — 13%, 
25%, and 40%. 




Fig. 10. Twenty-three neutrally buoyant ■ ellipses (p s = 1) at f a = 13%. {left) The velocity distribution, and (right) the 
angular distribution function. 


multi-ellipses: 

m= 

i = 1 

where N is the total number of particles and i indicates the *-th particle. The results in Fig. 10 show that 
the neutrally buoyant elliptic particles have a much higher probability to align in the flow direction than in 
the cross flow direction. 

To investigate the effect of the total particle volume fraction f a on migration, the simulations of 44 
and 70 elliptic particles are conducted, corresponding to f a = 25% and 40%; and Re& = 71.81 and 75.78, 
respectively. The results show that the maxima of the average particle volume fraction distribution across 
the channel become broader and the minimum gradually disappears as the total particle volume fraction f a 
increases (Fig. 9). It is observed that the rotational motion of an individual particle is also suppressed by 
the interactions from its neighboring particles. Thus, the lift due to rotation is reduced and particles may 
move more closely to the centerline. The velocity distribution of flow becomes flatter (cf. Fig. 11). These 
results suggest that the Segre-Silberberg effect becomes weaker due to multi-particle interaction, as the total 






Fig. 11. The velocity distributions of neutrally buoyant multi-particles at f a — 25% and f a — 40%. The velocity is 
normalized by the maximum velocity. 



Fig. 12. The averaged volume fraction distribution across the channel for elliptic particles with the densities of p# — 0.993 
and Rei = 109.38. and of p# ~ 0.9995 oA Rei, ~ 129.7. The total particle volume fraction for the two systems is f a 25% 


particle volume fraction increases. 


3.4. Non-neutrally buoyant multi-particles. When the particle density is slightly lighter than the 
fluid density (e.g., p s — 0.9995), the two maxima between the walls and the centerline on the averaged 
particle volume fraction distribution across the channel become much wider, as shown in Fig. 12, compared 
to that of the neutrally buoyant particles (Fig. 9). As the particle density further decreases to p, = 0.993, the 
two maxima finally merge together, thus the volume fraction distribution across the channel now has only 
one wide maximum around the centerline as shown in Fig. 12. The two shoulders on the curve of velocity 
profile become more prominent and the velocity profile becomes flat. The probability of ellipses aligning 
with the flow is also greatly reduced, as shown in Fig. 13. 

Finally, the system of 44 buoyant elliptic particles ( f a = 25%) with p s = 1.002 and p s = 1.015 are 
studied and compared with the system of 44 neutrally buoyant particles (p s = 1.0). When particles are 
slightly heavier than the fluid ( p s = 1.002), they move to the region closer to the walls than do the neutrally 
buoyant particles due to the inertial lift forces. In this case, the behavior of the multi-particle system is similar 
to that of a single particle. As shown in Fig. 15, the maxima of the averaged volume fraction distribution 
across the channel are located more closely to both walls for slightly heavier particles ( p„ = 1.002) than for 
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Fig. 13. The angular distribution functions for the same cases as in Fig . 12. 



Fig. 14. The velocity distributions for the same cases as in Fig. 12. 


neutrally buoyant particles. As the particle density increases to p, = 1.015, more particles move from the 
wall regions to the central region. By comparing the angular distribution functions of the neutrally buoyant 
particles and the heaviest particles (p s = 1.015) in Fig. 16, it can be seen that when p s = 1.015 the likelihood 
of particles aligning with the flow is greatly reduced by sedimentation, which turns ellipses from the flow 
direction to the cross flow direction. 

The inertial lift associated with the curvature effect in this case still exists, although it is relatively weak. 
Therefore, the averaged particle volume fraction distribution across the channel, the angular distribution 
function, and the velocity distribution, respectively shown in Fig. 15, 16, and 17, demonstrate the combined 
effects of the Poiseuille flow and sedimentation. It is clearly shown in Fig. 17 that the velocity profile for the 
heavier buoyant particles of p s = 1.015 is broader than that for the lighter buoyant particles of p„ = 1.002, 
and the velocity profile for the lighter buoyant particles of p s — 1.002 is broader than that for the neutrally 
buoyant particles. 

4. Conclusions. The lattice Boltzmann method has been employed to simulate the dynamic migration 
and rotational behavior of a single elliptic particle and multiple elliptic particles in the planar Poiseuille flow. 
The buoyancy effect on the lateral migration of the particles is analyzed. The total particle volume fractions 
of 13%, 25% and 40%, are considered in the simulations. Several conclusions can be drawn from this work. 

First of all, for the system consisting of a single neutrally buoyant elliptic particle or multiple particles. 
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Fig. 15. The averaged particle volume fraction distributions for elliptic particles with the density of p s ~ 1.0 at Re& ^ 
75.78, of p s — 1.002 at Re& — 82.81 and of p s — 1.015 at R e& — 114.06. The solid area fractions for the three systems are the 
same at f a — 25%. 



Fig. 16. The averaged angular distribution across the channel for the same cases as in Fig. 15. 


the Segre-Silberberg phenomenon is clearly observed in the simulations. A single elliptic particle migrates 
to an equilibrium position between the walls and the centerline of channel due to the inertial lift associated 
with the curvature of velocity profile. The equilibrium position of the particle in the cross flow direction 
is independent of the initial position and orientation of the particle. At the equilibrium, an ellipse is 
much more likely to align its long end along the flow direction than along the cross flow direction, because 
the torque from the flow on the particle is weaker when the particle align its long axis with the flow. 
Our simulations demonstrate that the Segre-Silberberg effect also exists in a multi-particle system, and is 
responsible for driving the particles away from the centerline. Consequently, the particle volume fraction 
distribution across the channel has two maxima located at both side of the centerline and has one minimum 
at the centerline. The velocity profile has two shoulders corresponding to the maxima in the particle volume 
fraction distribution. As the total particle volume fraction increases, the maxima in the particle volume 
fraction distribution become wider and the minimum disappears gradually. In the mean time, the shoulders 
in the velocity profile become more and more prominent and the velocity profile becomes more and more 
blunt. It is clear that the Segre-Silberberg phenomenon in a single elliptic particle system is caused by 
the wall lubrication repulsion, an inertial lift due to shear slip, a lift due to rotation and a lift associated 
with velocity curvature. However, the Segre-Silberberg effect is greatly reduced as the total particle volume 
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Fig. 17. The averaged velocity distribution across the channel for the same cases as in Fig. 15. The profile is normalized 
by the maximum velocity of suspensions for each case. 


fraction increases, and virtually vanishes as the total particle volume fraction increases beyond 40%, because 
multi-particle interaction changes velocity profile and reduces rotational motion of the particles. Ellipses 
preferentially align its long axis along the flow direction for the values of the total particle volume fraction 
considered here. 

Secondly, when the particle density is slight, ly heavier than the fluid, a single particle leads ahead the flow 
and moves to a position away from the channel centerline, because the inertial lift points to the wall. When 
the particle density is large enough, the particle moves back to the centerline. In this case, the slip velocity 
overshadows the effect of the velocity profile. The wall lubrication repulsive force becomes dominant and 
the particle motion is similar to sedimentation. The ellipse turns to the cross flow direction. The behavior 
of the two-way migration of the ellipse is the same as that found by Feng, Hu and Joseph [6, 7] and Huang, 
Feng and Joseph [10, 11] for a circular particle. This two-way migration behavior remains in a multi-particle 
system, i.e., the particles slightly lighter than the fluid move away from the centerline, and the particles 
which are heavy enough move back to the centerline. The likelihood that the ellipses align their long axes 
with the flow direction is greatly reduced due to the sedimentation effect,. 

And finally, when the particle density is smaller than the fluid density, a particle lags behind the velocity 
of flow, the inertial lift and wall repulsion grow and push the particle to the centerline. As the density dif- 
ference between the particle and the fluid (buoyancy) increases, the particle is concentrated at the centerline 
region. Oscillations about the centerline are observed in some cases. For multi-particle systems, the maxima 
in the particle volume fraction distribution across the channel increases, while both the angular distribution 
and the velocity profile become flat, as respectively shown in Figs. 14 and 1(3, due to buoyant forces and 
multi-particle interactions. 
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